Steady State of microemulsions in shear flow 
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Steady-state properties of microemulsions in shear flow 
are studied in the context of a Ginzburg-Landau free-energy 
approach. Explicit expressions are given for the structure 
factor and the time correlation function at the one loop level 
of approximation. Our results predict a four-peak pattern for 
the structure factor, implying the simultaneous presence of 
interfaces aligned with two different orientations. Due to the 
peculiar interface structure a non-monotonous relaxation of 
the time correlator is also found. 

PACS numbers: 61.20.Gy; 82.70.-y; 83. 50. Ax 



Self-assembling amphiphilic systems are binary or 
ternary mixtures with a surfactant forming interfaces be- 
tween the other fluid components. These systems show 
a very rich phase behavior with ordered and structured 
disordered phases at different temperatures and relative 
concentrations |l| . Particularly interesting for the appli- 
cations is the microemulsion phase where coherent do- 
mains of oil and water on scales between 100 and 1000 A 
form disordered isotropic intertwined structures on larger 
scales 

While a satisfactory comprehension of the equilibrium 
behavior of amphiphilic systems is now available [|l[ the 
dynamical features are far less understood. In particular, 
in microemulsions, the existence of correlated mesoscopic 
structures makes the effects of imposed flows very differ- 
ent than in simple fluids. Interfaces between oil and water 
are elongated in the flow direction and an excess stress 
appears in the system. This can give a peculiar dynamics 
with unusual relaxation properties and a non-Newtonian 
rheological behavior Moreover, from a statistical me- 
chanics point of view, complex fluids in external flow are 
an interesting example of an out of equilibrium driven 
system [Q. 

In this Letter we study the behavior of microemulsions 
in shear flow with a continuum free-energy approach sim- 
ilar to that used by Frcdrickson and Larson 1^ for the 
study of rheology of block copolymers and by Onuki and 
Kawasaki Q and others for the critical properties. 
A Ginzburg-Landau model for describing the rheological 
behavior of ternary mixtures was considered by Patzold 
and Dawson [|| . However the results of [D are based on 
an evaluation of the structure factor obtained numeri- 



cally or analytically only in the limit of vanishing shear 
rate 7. Here we solve exactly the model for any value of 
7 showing that in an intermediate shear range the struc- 
ture factor is characterized by four pronounced peaks in 
the plane of the shear and flow directions. This new phe- 
nomenon can be interpreted in terms of a complex spatial 
pattern where interfaces with two different orientations 
coexist; their relative abundance is tuned by 7. We also 
consider the behavior of the two-time correlator. Due to 
the coupling between the flow and the interface structure, 
this function has a remarkable non monotonic behavior 
characterized by the presence of maxima preceding a fast 
asymptotic decay. 

The Ginzburg-Landau free-energy generally used to 
describe the equilibrium properties of ternary mixtures 
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where the scalar field represents the concentration dif- 
ference between oil and water components. The amount 
of surfactant present in the system is related to the value 
of (70 which can be negative, favoring the appearing of 
interfaces. The term proportional to c makes stable the 
free energy at large momenta and weights the curvature 
of interfaces. The quadratic part of the free-energy in the 
disordered phase with a > and < Aac gives a suitable 
description of the microemulsion phase; in particular the 
space correlation function 
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well fits the experimental data |10|-|13|]. The characteris- 
tic length d = 27r(l/2-\/a/c — (7o/4c)^^/^ represents the 
size of coherent oil or water regions which are correlated 
up to a distance ^ — (l/2-y/a/c+ (jio/4c)~-^/^. At negative 
50, the structure factor obtained by Fourier transforming 
Eq. (2) has a peak at a finite value of the momentum 
given by /cm = -^/[goL/^c. Finally, the non-linear terms in 
eq.(l) describe the possible effects due to mode coupling 
and become more important in approaching the phase 
boundary. 



1 




The kinetic behavior of the mixture, neglecting hydro- 
dynamical effects 1 14], is described by the convection- 
diffusion equation [15 
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where the velocity field is given by 
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7 being the shear rate and Cx the unit vector in the flow 
direction. The stochastic term 77 is a gaussian white noise, 
representing thermal fluctuations, with zero mean and 
correlation {■q{r,t)'q{r',t')) = -2TTV'^5{r - f)S{t - t'), 
as required by the fluctuation-dissipation theorem, where 
(...) means the ensemble average. Here F is a mobility 
coefficient and T the temperature of the heat bath in 
contact with the fluid. 

The cubic terms appearing in the functional derivative 
of Eq. (3) will be treated in the one-loop approximation 
which describes appropriately disordered phases. The 
resulting equation for the Fourier components (p{k, t) is 
given by 
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where X{k) = cfc** -I- g^k^ +a, So{t) and S2{t) are defined 
self-consistently by the relations 
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and A is a phenomenological ultraviolet cut-off. Starting 
from Eq.(5) a straightforward calculation give the 

dynamical equations for the various correlation functions 
prf . We first focus on the properties of the structure fac- 
tor {(p{k, t)ip{—k, t)). Using the method of characteristics 
the t ^ 00 stationary expression of the structure factor 
C{k) can be calculated as 

C{k) = /'°"d,2FTi?2(z)e-r'^^'^^'(^-)[""+^"^'(^)+^^'(^)l 
"'0 

(8) 



where K{s) = {kx,ky + jk^s^kz), = 
a + b/2S^ + g^Sf, = go + .92^0° , 
are the limit for t ^ 00 of So{t) and 6*2 (i) 
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FIG. 1. The projections of the structure factor on the 
planes kz = (left column) and ky — (right column) at 
shear rates 7 — 0.25,2,4,8 (from top to bottom). The axis 
kx points towards the right. The vertical scale is the same in 
all the pictures. 

Since ^ and d are the physical quantities to be com- 
pared with experimental data, the model parameters are 
tuned Q with the following procedure: We first calcu- 
late the one-loop renormalization of the bare coefficients 
a and g without shear and express the correlation lengths 
^ and d in terms of the renormalized parameters. Then, 
for a given choice of ^ and d, we invert the above men- 
tioned relations and compute the values of a, go to be 
used in the case with shear. In the following we show re- 
sults for the case ^ = 7,d ~ 5. The other parameters are 
set to the values c^2,b = g2 = T= l. Similar results 
have been obtained for different choices of parameters 
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with d/^ in the range between 0.5 and 2.5. 

The effects of the flow on C{k) are shown in Fig. 1, 
where the projections of the three-dimensional structure 
factor on the planes ky = and — are plotted at dif- 
ferent 7. At 7 = the structure factor is isotropic and its 
shape on each cartesian plane is that of a circular volcano 
with radius kM- When shear is applied, the projection on 
the kx = plane gives the same results of the 7 = case 
because the velocity field is along the x direction. On 
the kz = plane, instead, this pattern is progressively 
distorted as the shear rate is increased. For small 7 (see 
the case 7 = 0.25 in the left column of Fig. 1) the edge 
of the volcano becomes elliptical and slightly depressed 
along the kx = —ky direction so that two weak symmet- 
ric maxima located at 6 = aictan^ky / kx) — 7r/4, 37r/4 
are developed. This is confirmed by a small 7 expansion 
of Eq. (^). As 7 is increased each of these two maxima 
is splitted into two peaks, as shown in Fig. 1 at 7 = 2 
and 7 = 4. By increasing the shear rate, the maxima 
located at fc^; ~ become comparatively more important 
while the other peaks rotate clockwise and decrease their 
amplitude linearly in 7 until they disappear. Indeed, in 
the limit 7 — > 00, since terms proportional to powers of 
■fkx damp the exponential term on the r.h.s. of Eq. (|^), 
only the maxima of C{k) with fc^; = survive. Fig. 1 
shows that this is already observable for 7 = 8. Then, 
by letting kx — in Eq. we find the peaks located 
at ky = ikj^i. The description is completed with the 
projection of the structure factors on the plane ky = 
shown in the right column of Fig. 1. Here two peaks at 
kx — 0,kz — ikM, are observed to become sharper and 
sharper as the shear is increased. 

A peak of C{k) is generally interpreted as the signature 
of a characteristic length in the system which is propor- 
tional to the inverse of its position. In this case, since 
the system is not isotropic, to each maximum one asso- 
ciates three lengths, one for each space direction. Due 
to the symmetry fc — > — fc only the peaks not related 
by reflection around the origin can be considered. For 
large shear the existence of a single couple of maxima at 
kx = signals that interfaces are aligned along the flow. 
In the transverse directions the characteristic lengths are 
the same as without shear. For intermediate values of 
7 the additional peaks at {kx, ky,kz) reveal the presence 
of interfaces oriented with an angle a = arctan(— fca;/fcy) 
with respect to the flow, besides those aligned along the 
X direction. As 7 is increased these features are pro- 
gressively tilted in the direction of the flow while their 
relative abundance diminishes, as suggested by the be- 
havior of the maxima with kx ^ Q previously discussed. 
The existence of a four-fold peaked structure factor has 
been also reported in scattering experiments on segregat- 
ing mixtures (l^Jl^ . For small 7 it is interesting to note 
that, although the microemulsion is almost isotropic, the 
depression of C{k) along kx = —ky indicates a slight 



predominance of interfaces directed against the flow at 
a ~ 37r/4. 
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FIG. 2. The stress tensor as a function of the shear rate 7. 

Stretching of domains requires work against surface 
tension and results in an increase A?; of the viscosity [po[ . 
The excess viscosity is generally defined as Ary = (Jxyh 
with the shear stress given by axy ~ — J j^^kxky{g^ + 

2ck'^)C{k) ||]. In fig. 2 axy is plotted as a function of 
7. Shear thinning is observed. For small 7, Uxy grows 
linearly with 7 as also found in ||^. This behaviour is 
observed in correspondence of a structure factor similar 
to that shown in Fig. 1 at 7 — 0.25. When C{k) devel- 
ops four maxima at kz — Q, (^xy keeps increasing with a 
much smaller effective exponent consistent with the value 
0.13. For very large 7, when the peaks at fc^; = alone 
survive, Uxy decreases to zero because the excess stress 
vanishes when the interfaces are aligned with the flow. 
A similar behavior is shown by the first normal stress 
Ni^SjS^ikl- kDig'^ + 2ce)C{k) il. 

Next we consider the dynamical properties of the 
microemulsion phase. In the steady state the two- 
time correlation function T>{k,k' ,t) — {(p(k,t)cp{k' ,0)) — 

D{k,t)5 {K{t) + k'^ lH] can be computed through 

D{k,t) ^ C ^-!:^ii'is)la-+y-KHs)+cK\s)]as 

(9) 

D is the product of two terms: the structure factor 
evaluated at the translated momentum K{t) and an ex- 
ponential term. For a given k — K{Q), K{t) moves along 
the ky direction upward or downward depending on the 
sign of kx, as represented schematically in the lower in- 
set of Fig. 3. Therefore, when K{t) crosses the edge of 

the structure factor, C (^K{t)^ shows a pronounced max- 
imum. This feature is peculiar to microemulsions where 
the presence of interfaces gives the patterns of structure 
factor discussed above and cannot be observed in simple 
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fluids. The exponential term in Eq. (g) contains a poly- 
nomial of t and behaves as exp(— r7^fc^ct^/7) for long 
times. However, due to the presence of negative coeffi- 
cients, the decay can be much slower at earlier times, as 
it can be observed in the upper inset of Fig. 3. In partic- 
ular one can show that an inflection point is developed at 
tan^ = kx/ky = —jt in the sectors with kxky < 0. Other 
inflection points come out, even in the other sectors, due 
to the negative value of g^. Such a slow decay can pre- 
serve the observation of the maximum of C (^K{t)j in the 
full time correlator (^. In Fig. 3 the behaviour of D is 
shown for three typical values of k. In case A, K{t) does 
not intersect the volcano of the structure factor so that D 
decays monotonously. When the edge of C{k) is crossed 
once or twice, as in cases B and C respectively, D is char- 
acterized by a corresponding number of peaks. This rich 
relaxation behavior has never been described before to 
our knowledge and could be searched for in experiments 
0. 
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FIG. 3. The behavior of the time correlation function of 
Eq. (9) is shown for the three values of the wave- vector corre- 
sponding to the positions A, B, C in the lower inset of the fig- 
ure. The elliptical shape in this inset represents schematically 
the contour plot of the structure factor in the plane kz = 0. 
Solid, dotted, and dashed lines correspond respectively to the 
cases C, B, and A. In the upper inset the exponential factor 
in the r.h.s. of Eq. (9) is plotted for the same three values ok 
k. 

In conclusion, we have presented explicit expressions 
for the steady-state structure factor and the time corre- 
lation function of a Ginzburg-Landau model describing 
mixtures of oil, water and surfactant in the microemul- 
sion phase under the action of a shear flow. The struc- 
ture factor exhibits a four peak pattern which implies a 
non-uniform orientation of interfaces with two preferred 
directions. The time correlator, reminiscent of the inter- 
face structure, shows maxima superimposed on a global 



decay. All these predictions can be tested experimentally 
and show that together with the equilibrium behavior 
also the dynamics of complex fluid is characterized by a 
very rich phenomenology. 

We thank Kenneth Dawson for helpful discussions. F.C. 
is grateful to M.Cirillo and R. Del Sole for hospitality 
in the University of Rome. F.C. acknowledges support 
by the TMR network contract ERBFMRXCT980183 and 
by MURST(PRIN 97). G.G. acknowledges support by 
MURST (PRIN 97). 



[1] 



[2: 

[3: 
[4: 

[6; 

[7: 

[8 

[9 
[W 

[11 
[12 

[13 

[14 

[15 
[16 



For a review see, e.g., G. Gompper and M. Schick, 
Self-assembling amphiphilic systems in Phase Transitions 
and Critical Phenomena, edited by C. Domb and J.L. 
Lebowitz (Academic Press, New York, 1994), vol.16; 
K.A. Dawson in Structure and Dynamics of Strongly In- 
teracting Colloids and Supramolecular Aggregates in So- 
lution, edited by S.-H.Chen at al. (Kluwer Academic 
Publ. 1992). 

P.G. de Gennes and C. Taupin, J. Phys. Chem. 86, 2294 
(1982). 

See, e.g., R.G. Larson, The structure and Rheology of 
Complex fluids (Oxford University Press, New York, 
1999). 

B. Schmittmann and R.K.P. Zia, Statistical Mechanics of 
Driven Diffusive Systems in Phase Transitions and Crit- 
ical Phenomena, edited by C. Domb and J.L. Lebowitz 
(Academic Press, New York, 1995), vol.17. 
G.H. Fredrickson, J. Chem. Phys. 85, 5306 (1986); G.H. 
Fredrickson and R.G. Larson ibid. 86, 1553 (1987). 

A. Onuki and K. Kawasaki, Ann. Phys. (N.Y.) 121, 456 
(1979). A. Onuki, J. Chem. Phys. 87, 3692 (1987). 
M.E. Gates and S.T. Milner, Phys. Rev. Lett. 62, 1856 

(1989) . 

G. Patzold and K. Dawson, Phys. Rev. E 52, 6908 (1995); 
G. Patzold and K. Dawson, J. Ghem. Phys. 104, 5932 
(1996). 

G. Gompper and M. Schick, Phys. Rev. Lett. 65, 1116 

(1990) . 

M. Teubner and R. Strey, J. Ghem. Phys. 87, 3195 
(1987). 

B. Widom, J. Ghem. Phys. 90, 2437 (1989). 

G. Gompper and M. Schick, Phys. Rev. Lett. 62, 1647 
(1989). 

A. Gappi, P. Golangelo, G. Gonnella, and Amos Maritan, 
Nucl. Phys. B 370 659 (1992); P. Golangelo, G. Gonnella, 
and Amos Maritan, Phys. Rev. E 47 411 (1993). 
We are also neglecting the motion of surfactant assuming 
that its relaxation is fast compared to that of the order 
parameter (p (see [^). 

For a review, see A. Onuki, J. Phys.: Gondens. Matter 9 
6119 (1997). 

U. Marini Bettolo Marconi anf F. Gorberi, Europhys. 
Lett. 30 349 (1995). 



4 



[17] Details will be published elsewhere. 

[18] K. Migler, C. Liu, and D.J. Pine, Macromol. 29 1422 
(1996). 

[19] F. Corberi, G. Gonnella, and A. Lamura, Phys. Rev. 

Lett. 83 4057 (1999). 
[20] A. Onuki, Phys. Rev. A 35 5149 (1987). 
[21] A^i behaves at small 7 like 7^. 

[22] Notice the presence of the delta function 5 (^X{t) + fc'), 

as opposed to the usual S{k + k'), due to the distortion 
induced by the flow [^. 
[23] It is interesting to observe that an oscillating behavior 
of the time correlation function of microemulsions has 
been also predicted in the case without shear when hy- 
drodynamic effects are considered (G. Gompper and M. 
Hennes, Phys. Rev. Lett. 73 1114 1994). 



5 



